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ABSTRACT 

In symbolic computing, a major bottleneck is middle ex- 
pression swell. Symbolic geometric computing based on in- 
variant algebras can alleviate this difficulty. For example, 
the size of projective geometric computing based on bracket 
algebra can often be restrained to two terms, using final 
polynomials, area method, Cayley expansion, etc. This is 
the "binomial" feature of projective geometric computing in 
the language of bracket algebra. 

In this paper we report a stunning discovery in Euclidean 
geometric computing: the term preservation phenomenon. 
Input an expression in the language of Null Bracket Alge- 
bra (NBA), by the recipe we are to propose in this paper, 
the computing procedure can often be controlled to within 
the same number of terms as the input, through to the end. 
In particular, the conclusions of most Euclidean geometric 
theorems can be expressed by monomials in NBA, and the 
expression size in the proving procedure can often be con- 
trolled to within one term! Euclidean geometric computing 
can now be announced as having a "monomial" feature in 
the language of NBA. 

The recipe is composed of three parts: use long geomet- 
ric product to represent and compute multiplicatively, use 
"BREEFS" to control the expression size locally, and use 
Clifford factorization for term reduction and transition from 
algebra to geometry. 

By the time this paper is being written, the recipe has been 
tested by 70+ examples from [T], among which 30+ have 
monomial proofs. Among those outside the scope, the fa- 
mous Miquel's five-circle theorem [2], whose analytic proof 
is straightforward but very difficult symbolic computing, is 
discovered to have a 3-termed elegant proof with the recipe. 

ACM Computing Classification: 1. 1.1 [Symbolic and Al- 
gebraic Manipulation]: Expressions and Their Representa- 



tion; G.4 [Mathematical Software]: Efficiency. 
General Terms: Theory; algorithm. 

Keywords: Conformal geometric algebra, Null bracket al- 
gebra, Geometric invariance, Symbolic geometric comput- 
ing, Geometric theorem proving. 

1. INTRODUCTION 

Using geometric invariants in symbolic geometric comput- 
ing has been an active research subject in symbolic and 
algebraic computation. Apart from the benefit of better 
geometric interpretability when compared with coordinates 
[13] . geometric invariants have a salient feature of reduc- 
ing the size of symbolic manipulation. This is particularly 
valuable because a major difficulty in symbolic computing 
is middle expression swell. 

In projective incidence geometry, the proofs of many the- 
orems by the method of biquadratic final polynomials [3] 
can be so elegant that only bracket binomials occur in the 
whole procedure. The area method [2j also shares this fea- 
ture in its ratio-formed proofs of many theorems, i.e., the 
numerators and denominators are monomials of areas. The 
advantages of both methods are assimilated into the Cayley 
expansion theory developed in [11], by which this feature 
of projective incidence geometry is extended to projective 
conic geometry. 

In Euclidean geometry, with the introduction of inner prod- 
ucts, syzygies among basic invariants become much more 
complicated. Using distances for geometric theorem prov- 
ing was proposed in [4], and further developed in |12j by 
including the inner products of extensors. In [8] the covari- 
ant algebra of Euclidean geometry, the so-called conformal 
geometric algebra (CGA), and its invariant subalgebra, the 
so-called null bracket algebra, was shown to provide a nice 
algebraic setting for Euclidean geometric theorem proving. 
Unfortunately, except for some sporadic special cases, none 
of these methods has ever shown any binomial feature. As 
a benchmark problem, Miquel's 5-circle theorem [2j, whose 
analytic proof is straightforward but very difficult symbolic 
computing, was found a proof of 14 terms in 2001 [9]. The 
proof was full of pairwise term reductions based on compli- 
cated syzygies of advanced invariants. 



Under such background it comes as appalling as can be an 



observation that symbolic computing and theorem proving 
in Euclidean geometry have a "monomial" feature, or more 
generally, a term preservation feature in the language of 
NBA, if the recipe in this paper is used. 

By the time this paper is being written, 70+ Euclidean ge- 
ometric theorems in [T] have been tested, and 40+ have the 
term preservation feature. In particular, 30+ theorems have 
their conclusions represented by monomials in NBA, and are 
kept as monomials till the end of the proof. Concrete ex- 
amples include all the examples published in [8]. Using the 
method in that paper only one example preserves its num- 
ber of terms in the proof. Using our recipe in this paper 
ALL examples preserve their number of terms, thus it is 
impossible to find any analytic proof that is more elegant. 



Table 1: Comparison of proofs with [8] 



Example in [8] 


Conclusion 


Proof in [8] 


New proof 


No. 1 


1 term 


3 terms 


1 term 


No. 2 


2 terms 


2 terms 


2 terms 


No. 3 


1 term 


3 terms 


1 term 


No. 4 


2 terms 


3 terms 


2 terms 


No. 5 


1 term 


4 terms 


1 term 


No. 6 


1 term 


4 terms 


1 term 



By preserving the number of terms the computing burden 
is transmitted from addition to noncommutative multiplica- 
tion. One may shake head as to any possible simplification 
by algebraic manipulation of multiplication in place of addi- 
tion. Well, in an invariant symbolic system there are syzy- 
gies among basic elements. In manipulating such elements, 
multiplication preserves geometry while addition breaks it 
up. Symmetries in multiplication provide the most econom- 
ical way of avoiding or employing syzygies. It is easy to 
change multiplication to addition: just recall how coordi- 
nates are introduced. Generally it is very difficult to change 
addition to multiplication: just recall Cayley factorization 
[14] in projective geometry. 

Our recipe for symbolic computing in Euclidean geometry 
is: (1) employ multiplication, or more accurately, the geo- 
metric product in Geometric Algebra, from the representa- 
tion of geometric objects on, (2) preserve the multiplication 
through subsequent algebraic manipulation using the prin- 
ciple "BREEFS" [TT], and (3) replace addition by multipli- 
cation using Clifford factorization - the Euclidean version 
of Cayley factorization. In (1) we need to invent two new 
devices for the representation by multiplication, called nul- 
lifying operator and reduced meet product. In (2) we need to 
adapt the previous global invariant bracket-oriented princi- 
ple to a shift invariant neighborhood principle. In (3) we 
need a device to explore rational Clifford expansions sys- 
tematically - pseudodivision in NBA. These novelties will be 
introduced in Sections 4 and 5, with various illustrations. 

Geometric Algebra [5] is a version of Clifford algebra favor- 
ing the universal usage of its multiplication, the geometric 
product, instead of addition. Hestenes' vision of Geometric 
Algebra in place of the more commonly used Clifford algebra 
in matrix or hypercomplex numbers form, is fully justified 
by our theorem proving practice: replacing addition by mul- 
tiplication and prolonging the multiplication ( "long geomet- 
ric product'), are the simplest means of avoiding syzygies 



because the geometric product already incorporates various 
syzygies of inner products and determinants into its struc- 
tural symmetry. 

The Geometric Algebras developed for the conformal model 
of Euclidean geometry, CGA [6] and NBA [7], will be intro- 
duced in Section 2 from the implementation point of view. 
The geometry of long geometric product in NBA will be ex- 
plained in Section 3. In the end of this paper, a 3-termed 
analytic proof will be provided for the benchmark problem, 
Miquel's 5-circle theorem. 

2. CONFORMAL GEOMETRIC ALGEBRA 
AND NULL BRACKET ALGEBRA 

In [6], [7], [8] there have been detailed introductions of CGA 
and NBA. In this paper we concentrate on the case of 2D 
geometry only. We always use a boldfaced number or letter 
to denote a vector. 

In 4D Minkowski space we fix a null vector e. A null vector is 
a nonzero vector whose inner product with itself is zero. We 
call e the point at infinity of the Euclidean plane. Any null 
vector linear independent of e is a point in the plane. Two 
null vectors represent the same point if and only if they differ 
only by scale. This representation of the Euclidean plane is 
conformal but not isometric. To obtain an isometric model 
we only need to fix the inner product of any point with the 
point at infinity, e.g. to —1, as any two linear independent 
null vectors have nonzero inner product. 

To describe and analyze Euclidean geometry with the con- 
formal model we need a suitable algebraic language. The 
symbolic version of Clifford algebra in [5] is an ideal tool in 
that it prefers the usage of multiplication to addition. The 
multiplication, called geometric product, conglomerates all 
geometric relations within itself and is geometrically mean- 
ingful. The other two versions, the matrix version and the 
hypercomplex numbers version, emphasize the linear nature 
of the algebra, i.e., care more for addition than for multi- 
plication. In symbolic form, more addition leads to more 
algebra, and more multiplication preserves more geometry. 
This justifies the gist "geometric" in Hestenes' Geometric 
Algebra. 

Geometric Algebra is the unique algebra generated from an 
nD inner product space by an associative product, called 
geometric product, satisfying the generating relation that 
the geometric product of any vector with itself is the in- 
ner product. The geometric product is always denoted by 
juxtaposition. 

Geometric Algebra is graded, the grade ranges from to 
n. For an element A in this algebra, the i-graded part is 
denoted by (A)i. When i — and 1 it is the scalar and 
vector part respectively. Elements of grade i form a subspace 
of dimension G l n . In particular when i = n, the n-graded 
subspace is ID. Fix a nonzero element /„ in this space, then 
for any other n-graded element A n , the coordinate of A n 
with respect to the basis I n is the bracket of A n : 

[A n ] = ^ = A n I-' = I~ x A n . (2.1) 



The geometric product of an element A with I n 1 is called 
the dual of A: A~ = AI^ 1 . In particular A„ = [A n ]. 

The geometric product of two vectors is composed of two 
parts, the 0-graded part and the 2-graded part. They are 
respectively the inner product and outer product, denoted 
by dot and wedge. The outer product is just the exterior 
product in Grassmann's exterior algebra. 

12 + 21 = 2(12} =2(1-2), 2 (1 A 2) = 12 - 21. (2.2) 

For three vectors, 

2(123)i = 123 - 321, 2 (1 A 2 A 3) = 123 - 321. (2.3) 

Conformal geometric algebra (CGA) [6] is the Geometric 
Algebra established upon the conformal model. In CGA for 
2D geometry, we can pick out two scalar-valued functions 
generated by a sequence of 2k null vectors 1, 2, ■ ■ ■ , 2k: the 
0-graded part 



(12- ■■(2k)) = (12---(2k)) , 



(2.4) 

called the angular bracket, and the bracket of the 4-graded 
part for k > 1: 



[12...(2k)] = [(12...(2k)> 4 ], 



(2.5) 



called the square bracket. The two kinds of brackets generate 
a ring called null bracket algebra (NBA) [7], IS]- 

NBA should be implemented by realizing the following basic 
properties: 

(1) Multilinearity of geometric product, multilinearity and 
(anti-)symmetry of inner product and outer product, (|2.2|l 
and (|2.3[) from left to right. 



(2) Null symmetry: for null vector 1, 

1231 = -1321 = 1(2 A 3)1, 11 = 1-1 = 0. (2.6) 

(3) Shift and reversion symmetry: for vectors 1, 2, ■ ■ ■ , k, 



[12- --k] = -[2---kl] =-[kl---(k-l)], 

(12---k> = (2---kl>= (kl---(k-l)), 

[12---k] = [k(k-l)---21], 

(12 •••k) = (k(k-l)---21). 



(2.7) 



(4) Dual symmetry: for element A and r-graded element B r , 

(2.8) 



1(1 A 2 A 3)~ = -(1 A 2 A 3)~1, 
(A~) = [A], [A~] = -(A), 

A~B r = (~l) r (AB r )~, ~~ = -1. 



(5) Points 1,2,3 are collinear: [el23] = 0; points 1,2,3,4 
are cocircular: [1234] = 0. Lines 12, 1'2' are parallel: 
[el2el'2'] = 0; they are perpendicular: (el2el'2') = 0. 



(6) Contraction by Grassmann-Pliicker syzygy: for vectors 
1,2, 3, 4, 5,1', 2', 3', 

[1234][51'2'3'] - [1235][41'2'3'] + [1245] [3l'2'3'] 



= 1345 21'2'3' 2345 11'2'3' 



(2.9) 



(7) Contraction by inner-product bracket syzygy: for vectors 
1,2,3,4,5,1', 

1' • 1[2345] - 1' ■ 2[1345] + l' ■ 3[1245] 



= 1' • 4[1235] 1' • 5[1234] 



(2.10) 



(8) Null expansion: for null vector 1, 

121 = 2(1-2)1, 

1231 = 2(1-2)31-2(1-3)21. 



(2.11) 



(9) Trigonometric quartet expansion: for null vector 1, 

-[123415 • r] = (1234) [15 - r] 
+ [1234] (15 • • • r), 

-(123415 - r) = (1234) (15 • • • r) 
- [1234] [15 ---r]. 

(10) Trigonometric sextet expansion: for null vector 1, 
1 , 



(2.12) 



^[12314516- -r] = [123145](16 ■ ■ • r) 
+ (123145)[16 ■■■!■], 
i(12314516 - r) = (123145) (16 ■■■ r) 



(2.13) 



- [123145] [16 r]. 

(11) Rational sextet expansion: for null vectors 2,3,5,6, 



--[123456] [2356] = 2 ■ 3[1256] [3456] 
+5 ■ 6[1236][2345], 

i (123456) [2356] = 2 3 [1256] (3456) 



(2.14) 



-5 ■ 6(1236)[2345]. 

(12) Rational octet expansion: for null vectors 1, 2, 3, 5, 

i[12341256][1235] = 2 3[1256] [125134] 

-2 ■ 5[1234][123156], 
i(12341256)[1235] = 2 3(1256) [125134] 

2 5[1234](123156), 

i [13241526] [1325] = (1 ■ 5) (2 ■ 5) [1324] [1326] 



(13241526)[1325] 



+(1 • 3)(2 • 3)[1524][1526], 
(l-5)(2-5)[1324](1326) 



-(1 ■ 3) (2 • 3)[1524](1526). 

(2.15) 

(13) Reverse of the expansions from (|2.12|l to (12. 15f) : from 
right to left. They are basic Clifford factorizations. 



The derivation of the properties is easy. In Section [3] (|2.12|l 
and (|2.13[) are discussed, in Section [5] (12.14[) and (|2.15|l are 

analyzed. 

3. THE GEOMETRY OF LONG PRODUCT 

Prolonging the length of elements in inner products and 
brackets (i.e., determinants) is not only a device of sim- 
plifying symbolic computing, but an indispensable means of 
representing basic geometric relations. We take a look at a 
planar angle and its algebraic representation. 

An oriented angle is just a 2D rotation. Z123 can be repre- 
sented by three points: the vertex 2, a point 1 on the initial 
ray, and a point 3 on the terminal ray without requiring 
that 1,3 be equidistant from 2. Without resorting to in- 
equalities, any rational description of angle Z123 by points 
1, 2, 3 is accurate only up to kn. The equivalent classes of 
oriented planar angles modulo n are called full angles |2]- 



Let Z123, Zl'2'3' be two full angles. They are equal if and 
only if tan Z123 = tan Zl'2'3'. In the conformal model, the 
ratio of [el23] to (el23) is exactly tanZ123. So Z123 = 
Zl'2'3' if and only if 



[el23] _ [el'2'3'] 



(el23) (el'2'3') : 



(3.16) 



-[el23e3'2'l'] = [el23](el'2'3') - (el23) [el'2'3'] = 0. 

(3.17) 

Essentially, geometric product el23 represents full angle 
Z123, its sine and cosine are respectively the square and an- 
gular brackets, its reverse angle is e321. el2e34 represents 
full angle Z(12, 34). The sum of full angles e!23, el'2 '3' is 
their concatenation el23el'2'3'. These explain (|2,12[) and 
(|2.13p as expansions of the sines and cosines of angle sums. 



Using null expansion (|2.1ip and trigonometric expansions 
(I2.12[) and (|2.13(l . we easily obtain trigonometric explana- 
tions of all square and angular brackets. For example if e 
does not occur in the following vectors, then 



(1234) 
[1234] 

(123456) = 
[123456] = - 



(^12^230(34(^41 



cosZ(123,134) 



0(120(23(^341^4 



isinZ(123, 134), 



0(120(23 0(34^45 d^edei 



0(120(230(340(450(560(61 



cos(Z(123, 134) 
+Z(145,156)), 

sin(Z(123, 134) 
+Z(145,156)), 



and for the general case, 

(M2 • • • i2I+2) 
_ 0(iii2^i2'3 ' ' ' ^■2i + l i 2l + 2^ i 2l + 2'l 



cos(Z(iii 2 i3, 111314) 



+ Z(iii 4 i 5 , iiisie) H h Z(iii 2 ;i 2 i+i, iii2i+ii2!+2)); 

[ilia • ■ • i2i+2] 

0(i 1 i2 t (i2i3 ' ' ' < ^'2l + l i 2l+2 < ^hl+2 i l 



sin (Z(iii 2 i3, 111314) 



+ Z(iii 4 i 5 , iiisie) + 1 



Z(iiiaiiai+i, iiiai+i 121+2)). 



Here di2 is the Euclidean distance between points 1, 2; 123 
denotes the oriented circle through 1,2,3 sequentially, and 
Z(123, 134) is the full angle from the tangent direction of 
123 to that of 134 at any point of their intersection. 

If e occurs then the explanation is only slightly changed. 
For example, 

(el2345) = -di 2 d 2 30(34d45Cos(Z123 + Z345), 
[el2345] = -di 2 (i23o!34d45sin(Z123 + Z345), 

(3.18) 

and the general case follows similarly. The power of long 
geometric product comes from its intrinsic geometric nature. 

4. THE POWER OF LONG PRODUCT 



Below we present two novel devices in NBA, the nullifying 
operator and the reduced meet product. They function as 
the bridge between Grassmann-Cayley algebra and NBA, 
thus allowing to employ the full power of Cayley expansion 
theory [TT] within Euclidean geometry. 

Let 1 be a null vector in 4D Minkowski space. For vector 2, 

1 



2Vi(2) 



212 



(4.19) 



is the nullification of 2 with respect to 1. When 2 is null 
then iVi(2) = (1 ■ 2)2 represents the same point 2. When 2 
is not null, then Ni (2) represents the null vector other than 
1 in the plane spanned by 1,2 if the metric of the plane is 
Minkowski, or just 1 if the metric is degenerate. 

The reduced meet product of two elements 2 A 3 and 2' A 3' 
modulo vector 1 is 



(2 A 3) Vi (2' A3') 



[122'3']3 - [132'3']2 
[1233']2' [1232']3'. 



(4.20) 



The second equality is modulo 1, i.e., the two sides differ 
by A 1 for a scale A. This product is a reduced form of the 
classical meet product of elements 1 A 2 A 3 and 1 A 2' A 3': 

(1 A 2 A 3) V (1 A 2' A 3') = [122'3']1 A 3 - [132'3']1 A 2 
= [1233']1 A 2' - [1232']1 A 3' 

(4.21) 

Proposition 1. [Null duality] Let 1' be a null vector, then 

1'(1A2A3)~(1A2'A3')~1' = l'l{(2A3)Vi (2'A3')}~l'- 

(4.22) 

Proof. In Geometric Algebra we have the duality relation 
A~ A B~ = (AV B)~, (4.23) 

so 

l'(l A 2 A 3)~(1 A 2' A 3')~l' 
= 1'{(1A2A3)~A(1A2'A3')~}1' 
= 1'{(1 A2A3) V(1A2' A3')}~1' 
= l'l{(2A3)Vi(2'A3')}~l'. 

Example 1. [See [8], Example 5] If three circles having a 
point in common intersect pairwise at three collinear points, 
their common point are cocircular with their centers. 




Figure 1: Example 1. 

This is the most difficult example in [8]. We use the same 
geometric scenario: remove the collinearity constraint from 
the hypotheses, compute the conclusion expression to see if 
the removed constraint comes out as a factor. 

Same as in [5], during the computing all intermediate factors 
are saved in a set for later analysis. They are marked with 
under braces and are removed from subsequent steps. 



Free points: 0, 1, 2. 

Centers: 4 = center(012), 5 = center(013), 6 = center(023). 
Removed hypothesis: [el23] = 0. 
Conclusion expression: [0456]. 

The center of circle 123 is 

center(123) = N e ((l A2A3)~). (4.24) 

Substitute the expressions of the three circle centers 4, 5, 6 
into the conclusion, we get 

2 3 [0456] = [0(0 A 1 A 2)~e(0 A 1 A 2)~(0 A 1 A 3)~e 

(0 A 1 A 3)~(0 A 2 A 3)~e(0 A 2 A 3)~] 

= [0e(0MA2)~(0AlA3)~e(0AU3)~ 

(0 A 2 A 3)~e(0 A 2 A 3)~(0 A 1 A 2)~] 

= [OeO{(l A 2) V (1 A 3)}~e0{(l A 3)V 

(2 A 3)}~e(0 A 2 A 3)~(0 A 1 A 2)~] 

= -2 (e ■ 0)[0123] 2 [01e03e(0 A 2 A 3)(0 A 1 A 2)] 

= 2 2 [01e03e023012] 

= 2 2 [01e0e3032021] 

= 2 2 (e • 0)(0 ■ 1)(0 ■ 2)(0 ■ 3)[el23]. 
- . ' 

Explanation of the computing: 
Line 1: substitution. 

Line 2: change of order by symmetries (|2.8[1 and (12.71) . 
Line 3: apply (I4.22|) to the first two pairs of meet products. 
Line 4: expand meet products, make null expansion (|2.1ip . 
Line 5: apply (|2.3[) . Only one term is generated. 
Line 6: change of order between neighboring 0's by (|2.6|l . 
Line 7: null expansion. 

Now consider the most typical geometric construction: the 
intersection of two circles. Let 123 and 12'3' be two circles 
represented by circumpoints. Their point of intersection, 
denoted by 123 n 12'3', refers to the point of intersection 
other than 1, or in the tangent case, tangent point 1 itself. 
When 1 = e it is the intersection of lines 23, 2'3'; when 
2 = e it is the intersection of line 13 and circle 12'3'. 

In [7] the intersection was expressed as a linear combination 
of three vectors: either 1, 2, 3, or 1, 2', 3'. In this paper we 
propose the following representation: 

123 n 12'3' = JVi((2 A 3) Vi (2' A 3')). (4.25) 

It is easily verified that (|4.25[) equals the expressions in [7]. 

In (|4.25[) the reduced meet product has two ways of expan- 
sion, either by separating 2, 3 as in the first line of (|4.20|l . or 
by separating 2', 3'. If we compute the geometric product 

4(123ni2'3')5 = i4{(2A3)Vi(2'A3')}l{(2A3)Vi(2'A3')}5, 

(4.26) 

then previously we simply expanded the meet products in 
the same way and multiplied them with 1. Since the geo- 
metric product is associative, not only can we expand the 
meet products in different ways, but we can freely change 
the order of the four pairwise geometric products in (|4.26[) . 



Furthermore, if the intersection has more than one construc- 
tion, e.g. it is where three lines meet, we can even change 
the representations by different pairs of lines for each meet 
product. By prolonging the length of the geometric product, 
we gain a lot of freedom for the realization of "BREEFS" . 

5. BREEFS 

BREEFS - "Bracket-oriented Representation, Elimination 
and Expansion for Factored and Shortest result", was first 
proposed in [TT] to control the expression size in bracket al- 
gebra. In [8] it was extended to null inner-product bracket 
algebra, the subalgebra of NBA generated by angular brack- 
ets of length 2 and square brackets of length 4. In this sec- 
tion we extend it to long geometric products. The best way 
to explain this principle is through some working examples. 



Example 2. In the plane two circles intersect at points 
1,1'. Two secant lines through them intersect the circles at 
points 2,3 and 2', 3' respectively. Then 22'//33'. 




Figure 2: (a): Original theorem; (b): Two con- 
straints removed. 

In [8] one hypothesis of the theorem was removed and a ge- 
ometric factorization of the conclusion was obtained. In [TD] 
two hypotheses were removed and a geometric completion 
was discovered. The computing was 3-termed in [8] and 2- 
termed in [TD]. Both required Clifford factorization, [S] fur- 
ther used circular transform. Below we present a 1-termed 
computing without Clifford factorization nor transform. 

The configuration in Figure 2b is constructed as follows: 

Free points: 1,2,3, 1', 2'. 
Intersection: 3' = l'2' n 131'. 
Conclusion expression: [e22'e33']. 

[e22'e33'] = 2 _1 [e22'e3{(e A 2') V r (1 A 3)} 
l'{(eA2')V 1( (1A3)}] 
= 2 1 [el3l'][e3l'2'][e22'e31l'2'] (5.27) 

= 2(e-2')[2'2e31l']. 
BREEFS in (Qrgfl : 

(1) The first meet product has two neighbors in the long 
geometric product: 3 and l'. Expansion by separating 1,3 
would simply replace the meet product by 1. Separating 
e, 2' would produce two terms and a Clifford factorization 
(12.14(1 must be used to get a monomial result. 

(2) Similarly, the expansion of the second meet product has 
to separate e, 2', as e is a neighbor by shift symmetry. The 
last step uses null symmetry (|2.6[l and null expansion (|2.11|l : 

[e22'e31l'2'] = [2'e22'e31l'] = [2'e2'2e31l'] 

= 2(e-2')[2'2e31l']. 



The result is not what we expected. The two brackets [el23] 
and [12l'2'] representing the two removed hypotheses are 
not in the final result. They can be produced by rational 
expansion (|2.14|l : 

1. , 1 ■ 3[e232'lfl2l'2'l - 2 ■ 2'[el23][13l'2'] 

- e31l'2'2 = l - - ^ : l - - -. 

2 L J [1232'] 

Example 3. [See [S], Example 3] Let 1', 2', 3' be points on 
sides 23, 13, 12 of triangle 123 respectively. Then circles 
12'3', 1'23' and l'2'3 meet at a common point 4. 




Figure 3: Left: Example 3; Right: Example 4. 

In [8] a collinearity constraint was removed in order to ex- 
plore the dependency of the conclusion upon the constraint. 
Below we remove all the hypotheses (the three collinear- 
ity constraints) to explore the equivalence of the conclusion 
with the hypotheses. 

Conclusion: three circles 12'3', 1'23', l'2'3 concur, i.e., the 
intersection 4 = 12'3' n l'23' is on circle l'2'3: 

[1'2'34] 
= [l'(12'3' n 1'23')2'3] 

= 2" 1 [1'{(1 A 2') V 3 / (1' A 2)}3'{(1 A 2') V 3 / (l' A 2)}2'3] 
= 2~ 1 [ll'2'3'][2l'2'3'][l'23'l2'3]. 

(5.28) 

In conclusion expression [l'2'34], l',2',3 are antisymmet- 
ric. Because 3 is irrelevant to 4 = 12'3'nl'23', only when 4 
is between l', 2' can the BREEFS principle take effect. Then 
similar to Example 2, neighborhood consideration leads to 
unique monomial expansions of the meet products. 

Discarded hypotheses: three "circles" 1'23, 12'3, 123' con- 
cur (at "point" e). This can be represented by the incidence 
of the intersection 1'23 n 12'3 and circle 123'. Simply by 
interchanging the primes over the same letters, we get from 
(|5.28p the same effective part of the discarded hypotheses: 

[12'3l'23'] = -[1'23'12'3]. (5.29) 

(|5.29[) discloses the intrinsic equivalence between the con- 
clusion and the hypotheses. 

Example 4. [Miquel's 4-circle Theorem] Four circles inter- 
sect sequentially at pairs of points (1,5), (2,7), (3,6) and 
(4, 8). If 1, 2, 3, 4 are cocircular then so are 5, 6, 7, 8. 

This is a typical theorem whose analytic proof using coordi- 
nates is difficult although straightforward. In [7] a 5-termed 
NBA proof was found. Below we present a 1-termed proof. 

Free points: 1,2,3,4,5,6. 



Intersections: 7 = 125 n 236, 8 = 145 n 346. 

[5678] = -[5768] 

= -2~ 2 [5{(1 A 5) V 2 (3 A 6)}2{(1 A 5) V 2 (3 A 6)} 
6{(1 A 5) V 4 (3 A 6)}4{(1 A 5) V 4 (3 A 6)}] 

= 2" 2 [1256] [1456] [2356] [3456] [51236341] 
= -2 2 (1 ■ 5)(3 • 6) [1234]. 

(5.30) 

Similar to Example 3, neighbors of 7, 8 in bracket [5678] 
should be changed to 5, 6 to make the best of the BREEFS 
principle. The rest are two simple null expansions leading 
to the conclusion [1234]. If using the original [5678], then 
the two meet products 

{(1 A 5) V 2 (3 A 6)}{(1 A 5) V 4 (3 A 6)} 

are neighbors in the long geometric product. They share two 
pairs of points: (1,5) and (3,6), and should be expanded 
by separating the same pair. The result is 2-termed, and 
Clifford factorization (12.14[) must be used to return to one 
term. 

Summary of "BREEFS" in NBA: 

(1) To get factored and shortest result, choose suitable al- 
gebraic representations, do eliminations and expansions in 
long geometric products according to neighborhood consid- 
eration. Immediate neighbors have topmost priority, then 
next immediate neighbors, and so on. (2) If neighbors of 
a meet product can be altered, then rearrange neighbors 
by relevance to the meet product. (3) If an outer product 
is common to neighboring meet products, expand the meet 
products by splitting common outer product. 

Below we discuss Clifford factorizations (|2,14[l to (12.15[) which 
are inverse of rational Clifford expansions. They are the 
most important devices to reduce the number of terms of 
bracket polynomials. In particular, the first formula of (|2.14|) 
is the MOST frequently used Clifford factorization. Its proof 
is straightforward substitution of the Cramer's rule of 4 (or 
1) with respect to basis 2, 3, 5, 6: 

[2356]4 = [3456]2 + [2456]3+[2346]5 [2345]6, (5.31) 

into bracket [123456] and then null expansion (|2.11J1 . All 
other rational Clifford expansion formulas are easily proved 
in this way. However, the proof does not provide any expla- 
nation as to why from the left to right of (|2.14|l one should 
add the factor [2356] , and why the right side should contain 
any of the factors on the left. 

(|2.14[) from left to right can be understood from the aspect 
of pseudodivision. In Example 2, [123456] in (|2.14|l needs to 
be pseudodivided by [3456] and [1236] to get the quotients. 
If 4 is the leading vector variable, then since divisor [3456] is 
linear with respect to 4, we can substitute the expression of 
4 by 3, 5, 6 into the dividend to get a remainder. Similarly, 
if 1 is the leading vector variable and [1236] is the divisor, 
we can substitute the expression of 1 by 2,3,6 into the 
dividend to get another remainder. The unique common 
pseudocoefficient for both pseudodivisions is [2356]. Thus 
the quotients can be obtained by substituting (|5.3ip into the 
dividend. 



To understand (|2.14[) from right to left is much more diffi- 
cult. However, once a factor say [2356] is discovered from 
the right, e.g., by polynomial factorization in homogeneous 
coordinates, then the other factor / can be easily fixed. 
Since / is a multilinear function of its vector variables, a pro- 
cedure similar to multilinear Cayley factorization [14] leads 
to the result that / equals [123456] up to a constant factor. 

6. ELEGANCE OF ANALYTIC PROOF 

Example 5. [Miquel's 5-circle Theorem] Let there be a five 
star with vertices l', 2', 3', 4', 5', and armpit points 1, 2, 3, 4, 
5. The circumcircles of the five triangular wedges meet se- 
quentially at shoulder points l", 2", 3", 4", 5" respectively. 
Then the shoulder points are cocircular. 




Figure 4: Miquel's 5-circle Theorem. 

There is another saying is that this theorem is due to W. K. 
Clifford, the inventor of Clifford algebra. Its analytic proof 
is straightforward but complicated symbolic computing. In 
this section we present a 3-termed beautiful proof with no 
short of elegance than traditional synthetic one. 

We use the same construction of the configuration as in [9]: 
first the armpit points 1,2,3,4,5 as free points; then the 
vertices as intersections of lines: 

1' = 23 51, 2' = 12n34, 3' = 23 n 45, 
4' = 34H51, 5' = 45nl2; 

finally the shoulder points as intersections of circles: 

1" = 11'2 n 55'1, 2" = 22'3 n 11'2, 3" = 33'4 n 22'3, 
4" = 44'5 n 33'4, 5" = 55'1 n 44'5. 

By symmetry we only need to prove [l"2"3"4"[ = 0. 

[1"2"3"4"] 

= -2- 4 [l{(l' A 2) Vi (5 A 5')}{(2' A 3) V 2 (1 A l')} 
2{(2' A 3) V 2 (1 A l')}{(3' A 4) V 3 (2 A 2')} 
3{(3' A 4) V 3 (2 A 2')}{(4' A 5) V 4 (3 A 3')} 
4{(4' A 5) V 4 (3 A 3')}{(1' A 2) Vi (5 A 5')}]. 

(6.32) 

Consider the expansion of the first meet product. Its imme- 
diate neighbors provide no hint. A next immediate neighbor 
is 2, suggesting the separation of l', 2. Similarly, for the sec- 
ond meet product, its next immediate neighbor 1 suggests 
the expansion by separating 1,1'. The two expansions result 
in three terms from the first two meet products: 

1{(1' A 2) Vi (5 A 5')}{(2' A 3) V 2 (1 A l')}2 

= [151'5'][231'2']1212 + [1232'][15l'5']12l'2 

[1255'][231'2']11'12 



= 2 (-1 • 2[151'5'][231'2'] + 2 ■ l'[1232'][15l'5'] 
-1 ■ l'[1255'][23l'2'])12 

= 5i2l2. 

Suddenly it appears that the first two meet products are 
simply removed from (|6.32|l : 

[1"2"3"4"] 

= -2- 4 gi 2 [12{(2' A 3) V 2 (1 A l')}{(3' A 4) V 3 (2 A 2')} 

3{(3' A 4) V 3 (2 A 2')}{(4' A 5) V 4 (3 A 3')} 
4{(4' A 5) V 4 (3 A 3')}{(1' A 2) Vi (5 A 5')}]. 

(6.33) 

Although irrelevant to the proof, It is interesting to note 
that after eliminating l',2',5' from its expression, the in- 
termediate factor gi2 equals 

-2(l-2)(l-5)(2-3)[el23] 2 [el25] 3 [el45][e234][e345]Si 3524 , 

where 

Si 3524 = (e-l)(e-4)[e235]-(e-l)(e-5)[e234] 

+ (e ■ 2)(e ■ 5)[el34] 

is twice the signed area of pentagon 13524. 

Because of the symmetry in the geometric constructions, in 
(|6.33p the expansion of the two meet products between 2, 3 
results in an intermediate factor (723, and the expansion of 
the two meet products between 3, 4 results in an intermedi- 
ate factor t/34. Of course by direct computing we obtain the 
same result. Now (|6.33[) is changed into 

h = P23334[1234{(4' A 5) V 4 (3 A 3')}{(l' A 2) Vi (5 A 5')}]. 

(6.34) 

In (|6.34[) , immediate neighbors of the meet products suggest 
two ways of expanding them simultaneously, each resulting 
in three terms: either separate 4', 5 and 5,5', or separate 
3, 3' and l', 2. The latter expansion has three terms because 

[123432] = 2 (3 ■ 4) [1232] = 0. 

To choose between the two options, consider the next imme- 
diate neighbors. The first meet product has next immediate 
neighbors 3, 1, and the second has 2,4. They suggest the 
separation of 3,3' and 1', 2 respectively. The benefit is im- 
mediate null expansion of long brackets: 

h = -[1255'][453'4'][123431'] + [15l'5'] [3454'] [12343'2] 

+ [1255'] [3454'] [12343' l'j 

= -2(3-4)[123l'][1255'][453'4'] - 2 (1 • 2)[15l'5']x 
x [2343'] [3454'] + [1255'] [3454'] [12343' 1']. 

(6.35) 

The first round of elimination finishes here. In the second 
round, all remaining constrained points are eliminated: 



[1231'] 


1' 


2 


3[el23] [el25] [el35], 


[1255'] 


5' 


1 


2[el25][el45] [e245], 


[2343'] 


3' 


2 


3[e234][e245] [e345], 


[3454'] 


4' 


3 


4[el35][el45] [e345], 


[151'5'] 


l' 


1 


5[el23] [el55'][e235], 


[el55'] 


5' 


-2" 


1 [el2e45] [el25][el45], 


[453'4'] 


3' 


4 


5[e234] [e235][e454'], 



[e454'] = -2~ 1 [e34e51][el45][e345]. 

Then 

h = (1 • 2) (3 • 4)[el25][el35][el45] 2 [e245][e345] 

{2 ■ 3[el23][e234][e235](4 ■ 5[el25][e34e51] 
+ 1 ■ 5[e345][el2e45]) + [12343'l']}. 

(6.36) 

(|6.36|) contains a matching with (|2.15[) for Clifford factoriza- 
tion, 123456 in the first formula for e54312 here: 

4 ■ 5[el25][e34e51] + 1 ■ 5[e345][el2e45] 

= -i[e543e512][el45]. 

Then (|6.36p . and consequently the conclusion, is reduced to 

[12343'l'] = 2~ 1 (2-3)[el23][el45][e234][e235][e543e512]. 

(6.38) 

Below we prove (|6.38p . 
[12343' 1'] 

= 2~ 2 [1234{(2 A 3) V e (4 A 5)}e{(2 A 3) V e (4 A 5)} 
{(2 A 3) V e (1 A 5)}e{(2 A 3) V e (1 A 5)}]. 

By neighborhood consideration, the first meet product is 
expanded by separating 4, 5, the last meet product is ex- 
panded by separating 1, 5: 

[12343'1'] = 2~ 2 [el23][e234][12345e{(2 A 3)V e 

(4 A 5)}{(2 A 3) V e (1 A 5)}e5]. 

(6.39) 

The meet products in (|6.39|l form the double-line type in 
Cayley expansion theory [11] : 

e{(2 A 3) Ve (4 A 5)}{(2 A 3) V e (1 A 5)}e 
= e{{(2 A 3) Ve (4 A 5)} A {(2 A 3) V e (1 A 5)}}e 
= {(2 A3) V e (4 A 5) V e (1 A 5)}e(2 A3)e 
= [el45][e235]e23e. 
„™ (6.40) 
Then (]6.38[l is equivalent to 

[12345e23e5] = -2 (2 ■ 3)[e543e512]. (6.41) 

It must be pointed out that without using Cayley expansion 
we can simply expand the meet products in (|6.39[) by sepa- 
rating 2, 3 simultaneously. The result is 2-termed which can 
be contracted to one term using Grassmann-Plucker syzygy 
([2.9]) . The proof still remains 3-termed. 

Now Miquel's 5-circle theorem is equivalent to algebraic 
identity (]6.41[l . Use null symmetry and shift symmetry to 
rearrange the sequence so that two 3's are separated by three 
vectors, then use the trigonometric quartet expansion and 
factorization (|2.12l) . and null expansion (|2 . 11 p to get 

[12345e23e5] 
= [345e32e512] 

= -2 ((345e) [32e512] + [345e](32e512)) (6.42) 
= -2 2 (2 ■ 3)((e543) [e512] + [e543] (e512)) 
= -2(2-3)[e543e512]. 

This finishes the proof of the theorem. Summing up, with 
(Vi) denoting a reduced meet product with respect to i, the 
proving of the whole theorem proceeded as follows: 



[1"2"3"4"] — » [l(Vi)(V a )2(V a )(V 3 )3(V8)(V4)4(V4)(Vi)] 
— ► [12(V 2 )(V 3 )3(V 3 )(V 4 )4(V 4 )(Vi)] 
— ► [123(V 3 )(V 4 )4(V 4 )(Vi)] 
— ► [1234(V 4 )(Vi)] 
— > 3-termed h in (I6.35[) 
— > 3-termed h in (|6.36[) 

— > 2-termed (EMI: [12343'l'], [e534e512], 

succeeded by 

[12343'1'] — > [1234(V e )e(Ve)(Ve)e(Ve)] 
— > [12345e(Ve)(V e )e5] 

> [12345e23e5] 

► [e534e512]. 
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